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Black holes hold a tremendous discovery potential. In this paper the extent to which the Event 
Horizon Telescope and its next generation upgrade can resolve their structure is quantified. Black 
holes are characterized by a perfectly absorptive boundary, with a specific area determined by intrinsic 
parameters of the black hole. We use a general parametrization of spherically symmetric spacetimes 
describing deviations from this behavior, with parameters controlling the size of the central object 
and its interaction with light, in particular through a specular reflection coefficient [ and an intrinsic 
luminosity measured as a fraction 7 of that of the accretion disc. This enables us to study exotic 
compact objects and compare them with black holes in a model-independent manner. We determine 
the image features associated with the existence of a surface in the presence of a geometrically thin 
and optically thick accretion disc, identifying requirements for VLBI observations to be able to cast 
meaningful constraints on these parameters, in particular regarding angular resolution and image 
dynamic range. For face-on observations, constraints of order 7 < 10~*,P < 107! are possible with 


an enhanced EHT array, imposing strong constraints on the nature of the central object. 


I. INTRODUCTION 


The gravitational interaction is tested to exquisite pre- 
cision in the weak-field regime [1]. Gravitational-wave 
astronomy is now providing access also to strong-field 
regions, in violently dynamical regimes, and therefore to 
entirely new and exciting tests of Einstein’s theory [2-7]. 
Thus far, gravitational-wave data is entirely consistent 
with all predictions of General Relativity (GR) [8, 9]. 

It is, however, widely accepted that the description of 
gravity — as provided by the classical equations of GR — is 
incomplete. In particular, the theory breaks down in black 
hole (BH) interiors [10]. In addition, spacetime warping 
at BH horizons introduces puzzles in the semiclassical 
treatment of free fields [11]. The resolution of these and 
other issues seems to require an improved theory, and 
corresponding different spacetime geometries. The scale 
and nature at which the new effects play a role is unknown, 
and could range from “soft” horizon scale changes in the 
geometry, to “hard” effects doing away with horizons 
completely, at least in a classical sense [12-16]. 

It is up to observations to collect data pertaining to 
strong-field regions, in order to constrain the nature of 
dark compact objects [7]. Fortunately, BHs are also the 
simplest macroscopic objects in the cosmos, making the 
search for new physics associated with the absence of 
horizons particularly appealing [7, 17]. Due to their very 
nature, it is impossible to ever have observational proof 
that horizons do exist [7, 18], but one can certainly quan- 
tify the statement that there is, or not, structure close 
to the Schwarzschild radius inducing deviations from a 
perfectly absorptive behavior. In the past, observations 
of supermassive BHs, in particular Sagittarius A*, were 
used to qualitatively push any putative surface to Planck 
distances away from the horizon [19, 20]. The argument 


is based on thermodynamic equilibrium between the cen- 
tral object and its — visible — accretion disc. However, 
these arguments neglected strong lensing and conversion 
to other channels [7, 21, 22], and need to be revisited. 

In addition to gravitational-wave astronomy, opti- 
cal/infrared interferometry and mm-wavelength very large 
baseline interferometry (VLBI) are now able to scrutinize 
the region around BHs with unprecedented accuracy [23- 
25]. Of special interest to us here is the groundbreaking 
Event Horizon Telescope (EHT) observations relating to 
images of BHs. Most notably, the Event Horizon Tele- 
scope 1.3 mm wavelength images of the supermassive 
black holes in M87 [23] and Sgr A* [25] have revealed 
bright, ring-like emission on scales of the Schwarzschild 
radius. Here, we wish to understand how these observa- 
tions can be used to determine intrinsic features of BHs, 
such as the size of their boundaries, or their perfectly 
absorptive nature. This will give us information about 
the experimental designs that maximize the capabilities 
of testing these foundational aspects. 

We use geometric units with the speed of light and 
Newton’s constant c=G= 1. 


Il. THE SETUP 


Obtaining images of alternatives to BHs is a necessary 
endeavour to test GR and the BH paradigm. A systematic 
study would ideally take into account several non-trivial 
aspects which, in practice, would involve including param- 
eters and functional degrees of freedom [21] describing (i) 
the internal (bulk and surface) properties of the central 
object, (ii) the spacetime around them, and (iii) any mat- 
ter propagating around and interacting with the central 
object. 


Analyzing this complex question step by step, starting 
with the simplest possible models and gradually adding 
additional features, has the advantage of allowing us to 
examine the different physical aspects involved in a con- 
trolled manner. Moreover, simple models are interesting 
on their own right. Typically, their simplicity implies that 
the associated deviations from general relativity are less 
subtle than in more sophisticated models. Hence, simple 
models can be excellent markers of the boundary between 
observable and unobservable for a given VLBI design. 


A. Parametrizing the central object 


The simplest possible situation that we can analyze con- 
cerns spherically symmetric geometries with no additional 
external matter backreacting on the spacetime. In the 
context of “hard” but localized changes to the geometry, 
we can further make the assumption that the external 
geometry is given by the Schwarzschild solution, 


d 2 
ds? = — fat? + a + rd. (1) 


In practice, these assumptions involve neglecting first the 
external and matter parameters (points ii and iii in the 
list above), focusing on the internal parameters of the 
central object. 

The assumption of spherical symmetry constrains the 
internal parameters of the central object: the size is 
uniquely determined by a single parameter, its radius R. 
We are mostly concerned with discerning structure at the 
horizon, and therefore looking for changes that occur on 
small scales. As such, we set the surface of the central 
object to be 


R=2M(1+0), (2) 


where we will be mostly interested in situations in which 
e <1. Any coefficient describing bulk and surface prop- 
erties must be a function of the radial coordinate only. 
Generating images involves the ray tracing of null rays. 
All the information that is needed to image a given cen- 
tral object boils down to characterizing the its interaction 
with light. This interaction can take place both at bulk 
and surface levels, although for light rays it is reasonable 
to assume that these interactions take place mostly at 
the surface level. We can ignore then for the moment 
the bulk parameters (which, in practice, is equivalent to 
assuming that the bulk is optically thick). In this simple 
situation, we just need three parameters: a radius for the 
central object R, a specular (or elastic) reflection coeffi- 
cient [, and an intrinsic brightness 6 describing a locally 
isotropic surface emission due to a non-zero temperature 
(the latter is proportional to the dimensionless ratio be- 
tween the re-emitted energy and incident energy, [ « B, 
used in [21]). This situation can be implemented by intro- 
ducing a boundary in the Schwarzschild spacetime, thus 
creating a new spacetime with parameters {M,¢,T, B}, 


and modifying the propagation of null rays according to 
these coefficients. 

In practice, we will thus be analyzing images of space- 
times with parameters {M,«¢,T, 8}, and comparing these 
with BH images. BH images are associated with the 
subset {M,e <1, =0,B=0}.' As the nature of the 
central object can be very different when changing the val- 
ues of the parameters involved, it is reasonable to expect 
that these images can be quite different. Both specular 
reflection on and intrinsic brightness at the boundary 
of the central object are nonexistent in BH spacetimes, 
which implies the possible existence of new image features 
intrinsically associated with these processes. We will see 
that this intuition is accurate, although whether these 
differences can be measured with a specific experimental 
setup is a more subtle question. While the Event Hori- 
zon Telescope collaboration has partially analyzed this 
issue for specific situations [26], our analysis is exhaustive 
regarding the parameter space described above. 


B. Interaction between the central object and the 
accreting material 


As the disc material falls onto the compact objects, 
its velocity as measured by locally static observers is 
increasing. In fact, one can calculate the center-of-mass 
energy with which it collides with material at the (static) 
surface. Take then a proton or electron in the accreting 
material of mass mo colliding with another one at the 
surface. The center of mass energy is 


Eom = mov2,/1— Juv U4) Ula) (3) 


—_ mg 

yy VEG (4) 
where U(), Uy) are the 4-velocities of the particles [27]. 
In the last step we assume that one particle sits at the 
surface of the central object, while the other is described 
by a conserved energy parameter FE. For a particle that 
falls from the inner edge of the accretion disc, presumably 
at the innermost stable circular orbit, E = 22/3. 

We see that for « < 107!°, the center of mass energy 
for proton-proton collisions goes beyond the electroweak 
scale ~ 200 GeV. The collision by-product therefore con- 
sists of all possible Standard Model particles compatible 
with the scattering process. Since for small € all the by- 
products are strongly lensed and made to interact again 
with the central object, it is reasonable to expect that it 
thermalizes, reaching equilibrium with all the Standard 
Model species: our central object behaves essentially as a 
Hawking-radiating BH [28-31], although at a temperature 


1 A BHis strictly associated with « = 0 but, as situations with e < 1 
contain deviations proportional to « from a BH, the ubiquitous 
presence of numerical errors makes these indistinguishable in 
practice. 


dictated by the feeding disc. For small enough e, the CM 
collisions occur at high energy. For objects with effective 
temperatures down to of order 100 GeV, all Standard 
Model particles may be treated as essentially massless 
whereas for temperatures smaller than this there will 
be phase space suppression for the heavy gauge bosons 
and top quarks. In conclusion, we argue that the ultra- 
compact object should be re-emitting a nearly isotropic 
radiation, which is composed of ~ 1% of photons, the rest 
are neutrinos, electrons etc. 

The previous argument strongly suggests that ultra- 
compact central objects will be dimmer than previously 
thought in the electromagnetic band (and brighter in the 
neutrino band for example). A more rigorous analysis 
needs to be done to decide on the fate of heavier species 
(which may be unable to exit the central object due to 
the large gravitational potential). In fact, for cold (as 
measured at infinity) central objects with kT, < mc’, 
an elementary particle of mass m is unable to reach large 
distances and must fall back. In fact, for Sgr A*, EHT 
infrared observations and the assumption that the object 
emits blackbody radiation indicate that the central object 
is colder than < 10* K [26]. In this case, only photons 
and neutrinos reach infinity and therefore the fraction of 
accreting energy coming out as photons increases from 
1% to a sizeable fraction of the incoming flux [32].? 

Recent work on the appearance of horizonless geome- 
tries assumed that the surface of the central object emits 
blackbody radiation at the same temperature as the inner 
accretion flow 8 x 101° K [33]. This assumption is, in 
our opinion, ad hoc and lacks a sound theoretical basis. 
Indeed, were the central object + accretion disc system 
in thermal equilibrium, the locally measured (Tolman) 
temperature should be [34, 35] 


(5) 


expressing the fact that heat is also subjected to gravity. 
This law would predict a hotter central object. The ensu- 
ing radiation, after redshifts caused by the gravitational 
potential, would give energy conservation: what comes 
in from the disc, comes out as radiation. In other words, 
as is clear from Eq. (7) in Ref. [33], the more compact 
the central object, the colder it is. Thus it is far from 
equilibrium with the accretion disc. 


C. Ray-tracing and accretion flow model 


Calculating images of the accretion disc and reflecting 
surface requires determining the trajectories of photons 
(null geodesics) within the background spacetime geome- 
try. In this work we use the general-relativistic radiative 


? We thank Ramesh Narayan for bringing the Sgr A* constraints 
to our attention [26]. 


transport code BHOSS [36-38]. The geodesic equations 
of motion are formulated as a system of eight coupled 
first-order ordinary differential equations 

Ca Sore (6) 
where x2" denotes the position 4-vector, k“ the (contravari- 
ant) photon 4-momentum, If, the Christoffel symbols, 
and an overdot denotes differentiation with respect to the 
affine parameter. Equations (6) are integrated numerically 
in BHOSS using an adaptive fourth order Runge-Kutta- 
Fehlberg method (with fifth order error estimate). The 
integration tolerance is generally set to ~ 1071'°, suffi- 
cient for the purposes of this work. Details regarding the 
geodesic initialisation may be found in [37]. 

We consider the source of illumination to be a 
geometrically-thin and optically thick equatorial accre- 
tion disc. The disc is formally infinite in extent, with an 
inner edge specified by the radius of the innermost stable 
circular orbit (ISCO). The accretion disc material motion 
is assumed Keplerian, with four velocity, u“, given by [39] 


u4 = ul (1, 0, 0, Qx) , (7) 
where u’ and the Keplerian angular velocity, Nx, are 


; we VM : 

Y= faa (8) 
Once (a, k#) are calculated for a photon, the energy 
correction factor is calculated as g = —k,,u". The station- 
ary and axisymmetric nature of the spacetime geometry 
gives rise to two Killing vectors, from which one obtains 
k, = —E and kg = L,, where E and L, define the energy 
and angular momentum (along the spin axis) of the pho- 
ton, respectively. The energy correction factor may then 
be written as g = u! (E — OxL,). 

Finally, one must prescribe a local emissivity for the 
accretion disc. Since the disc is planar and opaque, we 
can neglect the effects of absorption and accumulation of 
intensity along the ray, thereby simplifying the problem of 
radiative transfer to determining the intersection of each 
ray with the disc surface, wherein each ray is prescribed a 
local monochromatic intensity, 7, which is weighted by g. 
In this work we adopt the profile j(r) « r~” with n = 3, 
where r is the disc’s radial co-ordinate and n has been 
chosen to ensure a finite total flux from the disc. 


III. IMAGE FEATURES ASSOCIATED WITH 
SPECULAR REFLECTION 


Given that the generation of images is a linear process, 
we can analyze the consequences of having a reflecting 
surface, or one with an intrinsic brightness independently. 

Hence, for the moment the only relevant physical pro- 
cess will be reflection, and we will focus on the case [ = 1 
that maximizes the associated features. Lower values of T 
would just make these features dimmer, as we will discuss 
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FIG. 1: Images of spacetimes with and accretion disc 
and a spherical mirror in Schwarzschild (left) and flat 
(right) geometries, with i = 0° and e« = 5,1, 0.1, in that 
order, from top to bottom. Note that the brightness of 
the accretion disc is noticeably different for the 
Schwarzschild and flat spacetimes, as the emission profile 
of the disc adjusts to the background gravitational field. 
Also, images for flat spacetimes do not have an 
associated mass scale, but we are still using the same 
field of view to simplify the comparison between both 
situations. Note also the distinctive features of curved 
spacetimes, which give rise to multiple rings, whereas a 
mirror in flat space only provides a single ring, from 
specular reflection. 


later. We have generated a series of images for different 
values of the radius R, parametrized in terms of a dimen- 
sionless parameter ¢€ as defined in Eq. (2). Our results 
for the images of horizonless objects when illuminated 
by the accretion flow described above are summarized 
in Figs. 1-2. The observer is located at an angle 7 with 
respect to the disc. For i = 0° the observer is at the poles, 
seeing the disc head-on. For 7 = 90° the observer is on 
the disc plane. 

These images display rich features associated with dif- 
ferent trajectories of light. However, not all these features 
are present for all values of €. In fact, images become in- 
creasingly complex as € decreases. It is thus useful to start 
discussing the features in the simpler images obtained 
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FIG. 2: Images for spherical mirrors in Schwarzschild 
(left) and flat (right) spacetimes, with 7 = 60° and 
€ = 5,1, 0.1, in that order, from top to bottom. 


for large values of €, and then discuss the new features 
that appear successively as the value of this dimensionless 
parameter decreases. The details of these images can be 
understood with the help of Fig. 3, which traces the origin 
of geodesics reaching the observer, and can be grouped 
as follows. 


A. Specular ring 


Specular reflection at r = R results in a “specular” ring, 
associated with the specular (or elastic) reflection in the 
surface of the photons emitted from the accretion disc. 
This feature is in fact present for all values of €, as long 
as there is an accretion disc around the central object. 

The largest value that we have considered for the images 
is « = 5, for which R = 12M is twice the radius of 
the innermost stable circular orbit (ISCO). In this case, 
the obtained images are a superposition of the direct 
and reflected images of the accretion disc. Moreover, 
due to the large value of ¢, differences between images 
obtained for the Schwarzschild and Minkowski spacetimes 
are negligible in this case (see Figs. 1-2). 

Figure 1 also shows that, as € decreases, the direct and 
reflected images of the accretion disc become separated, 
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FIG. 3: Origin of image geodesics for Schwarzschild reflecting surfaces for a face-on observer, i.e., i = 0° (upper left) 
and 7 = 60° (lower left). Pairs of critical geodesics which constitute the inner and outer boundaries of the “photon 
rings” and “reflected rings” in the upper left panel, as seen therein along the axis x = 0 and y > 0. The observer is 

located at z ~ 104 M. Reflected photons are denoted by solid lines and coloured by the quadrant of the surface they 
encounter. These photons comprise the reflected rings. Photons which are not reflected are denoted by dashed lines 
and coloured by the region of the accretion disc from which they are produced. The zoom-in view shows that the 

n = 3 photon ring (dashed yellow) and the m = 4 reflected ring (solid green) are in close proximity to one another but 

do not overlap. 


while also being possible to clearly see the effect that 
gravitational lensing has on the specular ring for the 
Schwarzschild spacetime. Note that the position of the 
specular ring changes noticeably with €, as expected from 
its interpretation. 

We are including plots of the Minkowski case in this 
section as a check of the sensibility of our results. The 
problem of reflection from a sphere is well studied in flat 
space [40, 41]. For an observer looking face on from large 
distances, points in the (thin) accretion disc at a distance 
r from the surface of the sphere will be seen to come from 
a ring of radius L with 


L= Rcoos(7/4+ arctan T), (9) 
1 

T=5 1 = VIF Be + VBy/-1 + 40 — VIF BE) (10) 
2 

r=—,. (11) 


For R = 2M and r = 6M for example, one finds 
£~1.591M. The end of the accretion disc (i.e., points 
at infinity) corresponds to points x ~ 0, and a ring at 
L = V2/2M ~ 1.414M, thus overall one gets a bright 
circular ring of thickness 0.176M. The results above can 
be extended to an arbitrary observation inclination an- 
gle [40, 41]. It is interesting that the projected ring has 
maximum thickness not for a head-on observer (@ = 0), 
but for an inclination 6 ~ 8.5 degrees. 


B. Lensing/photon ring 


The only role that gravitational lensing plays in the 
images discussed so far is the distortion (relative to the 
Minkowski case) of the specular reflection of the accretion 
disc. However, as € decreases, gravitational lensing be- 
comes large enough that photons from the accretion disc 


can circle the spherical mirror from below and then hit the 
observer without ever touching the surface of the central 
object, as shown clearly also in Fig. 3. This leads to a 
secondary image of the accretion disc. One can calculate 
the minimum radius beyond which this secondary ring 
appears. The behavior of null geodesics was studied in 
detail by Chandrasekhar [42] (see also [43]). The lensing 
deflection angle of a photon can be expressed in terms of 
elliptic integrals. In particular, if we define the perihelion 
distance P, then the lensing angle 0 of a photon in the 
Schwarzschild background is: ° 


Q = /(P-2M)(P+6M), (12a) 
».. Q-P+6M 
—.. (Q-P+2M 
Xoo = arcsin Q-Pt6M’ (12c) 


@:= ff [K(k?) — F(Xo0,k”)] —m. (12d) 


One can easily calculate the critical radius for which the 
bending angle is 7/2 or larger, corresponding to photons 
from distant regions of the accretion disc being seen by 
the face-on observer. The critical radius is R ~ 4.65958M 
(e ~ 1.32979). 

This analysis of null geodesics suggests that, for « < 
1.33, there should be an additional ring in the corre- 
sponding images, which we call “lensing/photon” ring 
(additional details on this name below). Note that the 
critical value € ~ 1.33 indicates the threshold below which 
a distant observer is able to receive at least one photon 
that has circled the spherical mirror. However, for the 
lensing/photon ring to be manifest in images, enough 
photons from the accretion disc must circle the spherical 
mirror. Hence, it is reasonable to expect that, depending 
on the resolution used to generate and analyze the images, 
as well as the emission profile of the accretion disc, the 
lensing/photon ring should start to be visible for a value 
of € slightly below this threshold. This is compatible with 
our images (Fig. 4). 

As shown in Fig. 4, we start seeing the lensing/photon 
ring for « = 1.25. For « = 1.25 the lensing/photon ring 
is quite dim. This changes if we consider smaller values 
of €, as the lensing/photon ring becomes brighter as the 
compactness increases, until this brightness hits a plateau 
for « < 0.5. However, the position of the lensing/photon 
ring does not change with ¢«. This suggests that this ring 
is associated with the geometry around and outside the 
spherical mirror, but not with reflection. To see this 
explicitly, we generated images with perfectly absorptive 
boundary conditions, verifying that such lensing ring is 
still present (see also Fig. 3). 


3 To keep in line with standard definitions — and in contrast with 
Chandrasekhar, we adopt the notation where the argument of 
the elliptic functions is k? instead of k. 
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FIG. 4: Cross-section of images for € = 1.25, « = 1.0 and 
€ = 1071, respectively (i = 0°). The lensing/photon ring 
starts to be visible for « = 1.25 (note the small peaks in 
between the accretion disc and the specular ring in the 
intensity profile), and it becomes brighter for decreasing 
values of «. As € decreases, there also appear additional 
rings between the lensing/photon and the specular rings. 


Regarding the name lensing/photon ring, we are par- 
tially following the notation in Ref. [44]. It is important 
to keep in mind that a general image (of a BH or any 
other object) does not have a single ring, but a sequence 
of them. What is generally called “photon” ring actually 
refers to a ring structure that shows different features 
associated with different photon orbits that can circle 
the central object multiple times. The authors of [44] 
proposed a specific way to split these features, further 
distinguishing between the lensing and photon rings, the 
first one being associated with photons that circle the 
central object once, thus covering a total angle of 27, 
while the second one is associated with photons that at 
least cover an angle of 57/2 and thus can circle the central 
object an unlimited amount of times. 

For our purposes, the most crucial aspect to keep in 
mind is that not all the features of the lensing/photon 
ring appear for the same values of compactness. In fact, 
in the terminology of [44], our results show that the fea- 
tures associated with the lensing ring appear for values 
of the radius greater than 3M (e = 0.5), in particular for 
RS 4.6 (€ S$ 1.3). We have discussed previously that, 
for this critical value of the radius, photons can start to 
circle around the spherical mirror and still reach a distant 
observer. As € decreases, photons can complete additional 
orbits around the spherical mirror before reaching a dis- 
tant observer, which means that additional photon orbits 
contribute to the lensing/photon ring. As a result, the 
lensing /photon ring gradually becomes as bright as the 
one in a BH image, and in fact becomes indistinguishable 
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FIG. 5: Profile that results from subtracting the 
cross-section for absorptive boundary conditions to the 
cross-section for reflective boundary conditions, for 
i= 0°, e= 107! and different values of [ = 1 (orange) 
and T = 0.3 (white). The features shown are associated 
with specular reflection and their size scales linearly with 
T, with the innermost ring corresponding to the specular 
reflection of the accretion disc. 


from the latter for « = 0.5. 
In conclusion, 


e The existence of a photon ring in an observational 
image is a relatively loose indicator of the compact- 
ness of the object, which can be anywhere in the 
interval e $1.3 (RS 4.6). 


¢ More refined evaluations of the intensity profile and 
substructure of the photon ring are needed to argue 
to further constrain the radius of the central object 
down to « $0.5 (R S 3.0). 


C. Reflected rings 


Aside from the specular and lensing/photon rings dis- 
cussed above, the spherical mirror leads to the appearance 
of additional rings associated with reflection, as can be 
seen for instance in Fig. 4. Recall that the specular ring 
is also associated with reflection, as it is the direct specu- 
lar image of the accretion disc. On the other hand, the 
lensing/photon ring is associated with photon orbits that 
circle the central object, without touching the latter (thus 
never being elastically reflected). The reflected rings are 
associated with photon orbits that both circle the cen- 
tral object but have also been elastically reflected in the 
latter.* 


4 Our numerical studies indicate that for « > 0.5 (i.e., outside the 
unstable photon orbit radius), the photon rings and reflected 
rings occur in pairs. Following Sec. IIIB, there are no photon 
rings (n = 0) and only the specularly reflected ring (hereafter 


To isolate the features associated with reflection, we can 
subtract the intensity profile for reflective and absorptive 
boundary conditions. Fig. 5 shows the result of this 
procedure for « = 1071. 


IV. IMAGE FEATURES ASSOCIATED WITH 
RE-EMISSION OR INTRINSIC BRIGHTNESS 


So far we have been analyzing situations in which the 
central object has a vanishing intrinsic brightness B = 0 
due to it being extremely cold. As we discussed above in 
Sec. IIB, the interaction between radiation and infalling 
material from the disc and the central object will heat the 
central object and lead to emission of all standard model 
particles. In addition, specular reflection (the channel we 
have studied so far) is only good as long as the wavelength 
of the radiation is larger than the typical particle separa- 
tion in the scattering surface. In other words, the surface 
of the reflecting object must be smooth on a wavelength- 
scale. Hence, including re-emission under the form of an 
intrinsic brightness is physically well-motivated. 

Fortunately, the linear nature of the problem simplifies 
the analysis. From the perspective of images, a nonzero 
value of B translates into the existence of a filled circle, 
with its brightness being a function of this parameter. 
It is possible to treat B as a completely independent 
parameter, to be constrained by observations. However, 
it is useful to relate B to the total flux radiated from the 
accretion disc, as described in the section below. 


A. Flux from the central object 


The accretion disc has a local emissivity profile 7 « r~”. 


We can formally calculate the total flux radiated from the 
accretion disc as: 


Frot x 2 [fio drdé = j(r) 4ar dr 
Disc "ISCO 


_ 40 (r1sco)?—” 


n—-2 , 


(13) 


where the factor of 2 accounts for emission from both the 
upper and lower parts of the disc. In this work we as- 
sume the background spacetime is Schwarzschild, yielding 
risco = 6M, and the index n = 3. The dimensionless 
total flux of the accretion disc in this work is therefore 
Fo, « 27/3. Since we omit absorption and are only con- 
cerned with emission at the disc surface, the constant of 
proportionality is arbitrary and hereafter fixed to unity. 


m = 0) when € = 1.3298. For « = 0.7603, the n = 1 photon 
ring and an m = 1 reflected ring are present, and this pattern 
continues as « + 0.5. The effect of decreasing epsilon is to reduce 
the diameter and thickness of all reflected rings, whilst leaving 
the geometry of the photon rings unaffected. See also Fig. 3. 


0.0 0.0 


FIG. 6: Images of spacetimes where specular reflection 
takes place, but with partial absorption (T = 0.5) and an 
intrinsic brightness (7 = 10~?) included. We take i = 0°, 
€ = 10~°, without filter (top row) and a Gaussian filter 
with the EHT angular resolution of 20 was (bottom row). 
We see that these values of 7 change appreciably the 
structure of the central depression in brightness. 


In considering the emission from the reflecting surface, 
we apportion a fraction, 7, of the total disc flux to the 
surface of the sphere, Fsphere = 7 Ftot = 27/3. From 
this, we uniformly distribute the aforementioned flux over 
the reflecting surface, yielding the local emissivity of the 
reflecting surface as: 


Jsurt = ch : (14) 
Note that Eq. (14) specifies that the local emissivity of 
the reflecting surface decreases (i.e., it becomes dimmer) 
as its radius increases, as expected. The reflecting surface 
is assumed to be stationary and hence the local energy 
correction factor is unity and the flux may simply be 
added to that from the accretion disc. 

We treat 7 as a phenomenological parameter, focusing 
on understanding the observability of image features for 
a range of its values. Specific proposals for the nature of 
the central object would result into specific values of 7, 
although more complete models of the central object, as 
well as its interaction with matter and light, are necessary 
to determine these values. Also, note that 7 relates the 
local emissivity of the surface and the total flux radiated 
by the accretion disc, and not the energy that the surface 
receives from the accretion disc. This should be taken 
into account when estimating the value of 7 from first 
principles (see the discussion below). Fig. 6 discusses the 


image features associated with re-emission. 

The parameter 7 encodes information about the physics 
of the central object, which receives direct radiation from 
the disc but also the accreting material itself. A simple 
estimate for thin discs shows that the accreting material 
contribution, M, is roughly one order of magnitude larger 
than the EM luminosity [45]. Hence, in our simplified 
model in which the accretion of matter is neglected, we 
are underestimating 7 by an order of magnitude. 

If a stationary state is reached, the central object should 
re-radiate all the incoming energy. However, as we dis- 
cussed in Section IIB, a substantial fraction may not 
emitted as EM radiation but in other channels °. If we 
consider that 1% is emitted as photons, then within our 
model one might expect 7 ~ 10-7. However, the effec- 
tive temperature as seen from infinity is expected to be 
lower than the disc’s, which could lower 7 by one order 
of magnitude. On the other hand, taking into account 
the accretion of matter would increase this by an order 
of magnitude, 7 > 1077. 

The above expectation holds provided there has been 
enough time to reach equilibrium, which depends on 
the value of «. As shown in Ref. [7], this happens for 
€Tsalpeter/(9.3M) > 1, with Teatpeter ~ 5 X 10’ yr the 
Salpeter time. Lowering the value of € can therefore com- 
promise reaching equilibrium in reasonable timescales. 
Moreover, the estimate for the timescale above assumes 
that no incident energy is spent in exciting the bulk de- 
grees of freedom of the central object [21, 22]. Hence, our 
discussion of the expected values of 7 above only apply to 
a subset of the available parameter space, but nevertheless 
provide a first estimation that can be used as reference in 
future studies, as well as combined with lower-bound con- 
straints that discard complementary regions in parameter 
space [50]. 


V. OBSERVABILITY OF IMAGE FEATURES 


In the previous sections, we have discussed the rich fea- 
tures associated with specular reflection and re-emission 
of incident energy. These images can be understood as 
the idealized limit in which observational errors vanish. 
The different sources of errors of VLBI observations limit 
our ability to detect such features in practice. 

The reconstruction of images according to a set of 
criteria, including systematic errors, but also aspects such 
as the array of telescopes used, is a complex problem in 
itself [51, 52]. However, it is standard to consider the 
application of a Gaussian filter as a first approximation, 
with the standard deviation o providing a measure of the 
angular resolution (see [53] for a detailed discussion). We 
will impose that the full width at half maximum, namely 


5 For example as positrons or high energy y rays, thus stimulating 
the calculation of accurate branching ratios, given the detection 
of anomalous rates from the galactic center [46-49]. 


2/2 1n 20, equals an angular resolution of 20 yas for EHT, 
while we consider tentative values of angular resolution of 
10 pas and 5 pas for ngEHT. The corresponding images 
are shown in Figs. 6 - 10. 

In general, we see that the image features associated 
both with specular reflection and re-emission are mostly 
filtered out for these values of angular resolution and 
that, generically, relatively large dynamic ranges (contrast 
between the brightest and dimmest points in a given 
image) are required to pick up these features. 

For an angular resolution of 20 yas and a dynamic 
range of ~ 10, face-on observations would allow us to 
constrain 7 < 107%, but do not constrain in practice the 
specular reflection coefficient [ (Figs. 6 and 7). Higher 
inclination angles translate into less stringent (or virtually 
nonexistent) constraints (Fig. 8). Lowering the angular 
resolution to 10 as does not change much the situation 
(Fig. 9). However, if the dynamic range increases to ~ 100, 
which is within the planned capability of a future ngEHT 
array [54], face-on observations would allow us to lower 
the constraint on the re-emission parameter by one order 
of magnitude, 7 < 10-4, and also constrain the specular 
reflection coefficient [ < 10~!. Again, higher inclination 
angles would weaken these constraints (Fig. 10), owing 
to much greater image flux asymmetry and the choice of 
azimuthal cut made to obtain the profile cut. 


VI. RECENT COMPLEMENTARY ANALYSES 


There has recently been progress in similar, complemen- 
tary directions. Naked singularities have been studied, 
from an imaging point of view, in Refs. [55] and [56]. 
The latter paper focuses on the contribution of higher- 
order lensed images. Being naked singularities, the fate 
of photons that reach the core is unclear; likewise, geome- 
tries with stable light rings are expected to be unstable, 
thus raising questions about the viability of such mod- 
els [57, 58]. Similar transparent matter is assumed to 
compose wormhole-type geometries, in a recent study of 
their shadow polarization pattern [59]. As we discussed, 
the interaction with the matter composing the central ob- 
ject is a fundamental ingredient in their appearance and 
imposes stringent constraints on the underlying physics. 

Other situations in which the interiors of the central 
objects are transparent have been analyzed in [60] and [61— 
63]. The image features obtained in these cases are asso- 
ciated with photons that circle the central object but also 
travel across it, and are thus complementary to the ones 
analyzed here, as we are assuming the opposite situation 
in which the central object is optically thick. This consid- 
eration translates into the energy-conservation constraint 
1=x«x+T+T, where x is the absorption coefficient, I’ 
the specular reflection coefficient, and [ « B is the di- 
mensionless ratio between the incident and re-emitted 
energy. While we find the assumption of optical thickness 
of the central object to be physically well motivated, we 
understand the value of adopting an agnostic perspective 
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FIG. 7: Results of applying a Gaussian filter with the 
EHT angular resolution of 20 pas for i = 0°, € = 1073 
and [ = 1, 7 = 0 (top row), fr = 0.5, 7 = 10-3 (middle 
row) and [ = 0, 7 = 0 (bottom row). We see that, even 
for angular resolutions that cannot resolve the internal 
structure, large enough dynamic ranges can allow us to 
distinguish between the situations in the third row and 
those in the first two rows. However, with this angular 
resolution it is generally difficult to disentangle the 
effects associated with the parameters [ and 7, as the 
similarities between the figures in the first two rows 
illustrate. 


and analyzing the observability of the features associated 
with (partially) transparent central objects. The image 
features analyzed in [60] and [61-63] would appear in situ- 
ations in which «+ T+T < 1, with a strength proportional 
to the magnitude of this deficit, and would be added to 
the ones analyzed here (which are present whenever I’ or 
T x B are non-vanishing). 


VII. CONCLUSIONS 


VLBI observations have the potential of unveiling the 
structure of BHs, and confirming whether general rela- 
tivity describes it properly. However, this requires that 
predictions of alternative models in which this structure 
is substantially modified, or even replaced altogether by 
something else, are falsified. In this paper we have used a 
previously proposed general parametrization of spacetimes 
that, besides the perfect absorption characteristic of BHs, 
can describe the specular reflection and re-emission of 
the incident energy, which can be characteristic of exotic 


ae = 1.0, 7=0.0 


ae = 0.0, 7=0.0 


5 0 5 W 15 20 


x (M) 
FIG. 8: Same as Fig. 7, but for 7 = 85°. We see that the 


difference between these situation shrinks, until 
disappearing in practice, as inclination increases. 


compact objects. We have presented a detailed analysis 
of the image features in the presence of a geometrically 
thin and optically thick accretion disc, and determined 
how experimental limitations impact the observability of 
these features. 

Our main results regarding image features are the fol- 
lowing: 


e Specular reflection manifests in an additional ring 
structure if the central object is compact enough. 


e Re-emission manifests in a central region with uni- 
form brightness. 


These two features are nonexisting in BH spacetimes, and 
thus provide a concrete and novel way of testing whether 
a given VLBI source is a BH. 

However, we have shown that, for the image dynamic 
range and angular resolution characteristic of EHT, and 
the ideal situation of 7 = 0°, it is possible to constraint 
the re-emission channel (7 < 1073), but not the specular 
reflection channel. We have also shown how improve- 
ments in image dynamic range and angular resolution 
such as the ones expected to be achievable in ngEHT can 
noticeably change the situation, leading to more stringent 
constraints re-emission channel (7 < 1074) and specular 
reflection channel (T < 10~*), at least for i = 0°. This 
provides further motivation for the improvement of VLBI 
observations, as well as theoretical modeling aimed at 
extracting more precise predictions. As we argued in Sec- 
tion IV A, known physics suggests that 1 > 10~? once the 
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FIG. 9: Results of applying a Gaussian filter with an 
angular resolution of 10 yas (top row) and 5 pas 
(bottom row) for i = 0°, «= 1073, f = 0.5 and = 107%. 
For 10 pas, the structure of image features associated 
with specular reflection cannot be discerned but, if the 
dynamic large is large enough, it is possible to constraint 
the surface parameters [ and 7. Only the most 
optimistic value of angular resolution (5 jas) can pick up 

the innermost structure of the simulated image. 


accretion of matter is accounted for, if € is greater than 
10~'° (M87) or 10~'8 (SgrA). Hence, ngEHT will con- 
strain significantly the nature of the dark central object 
even at such small scales. 
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features associated with the existence of a surface. 
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